On the Reception of Binary Signals in 
the Presence of a Small Random Delay* 
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(Manuscript received March 22, 1965) 

Receiver design for a binary communication system which operates over 
a linear channel with a random delay is considered. It is assumed that the 
variance of the random delay is very small and that the rate of growth of its 
moments is restricted. Under certain smoothness requirements on the re- 
ceived signal an approximation to the test statistic, which is optimum in the 
Neyman-Pearson sense, is derived for the case of gaussian receiver noise 
with covariance R(t) = R (0)e~ p ' T '. It is found that the test statistic, 
which in general is nonlinear, assumes the linear form of a crosscorrelator 
when phase reversal signaling is employed. 

The case where the noise is white and phase reversal signaling is used is 
investigated. The correlation waveform in this case is found to consist of the 
expected value of the received signal plus a term dependent on the slope of 
the signal when the delay is equal to its mean value. 

I. INTRODUCTION AND SUMMARY 

In any practical communication system the signal arrival time is 
never exactly known. This results in a degradation of the average sys- 
tem performance. It would be of considerable interest to determine the 
receiver which minimizes the effect of this uncertainty on system error 
performance. A special case of this problem will be considered here. 

Helstrom 1 has studied the detection of signals of unknown arrival 
time using the method of maximum likelihood with particular emphasis 
on the radar problem. Brown and Palermo 2 consider system performance 
in the presence of random delays with applications including least 
squares filtering and sampling with time jitter. Balakrishnan 3 and other 

* This work is based on a portion of a thesis entitled "Binary Signal and Re- 
ceiver Design for Linear Time Invariant Channels". This thesis was accepted by 
the faculty of the Graduate Division of the School of Engineering and Science of 
New York University in partial fulfillment of the requirements for the degree of 
Doctor of Engineering Science, Oct. 1964. 
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researchers have also considered the problem of time jitter in sampling. 
However, to the best of this authors knowledge, no optimum statistical 
test, or approximation thereof, has been determined for detection in 
the presence of a random delay. 

In the subsequent analysis we investigate binary communication for 
the case where the variance of the random delay is "very small." It is 
assumed that the transmitted signals and the channel impulse response 
are such that the received signal satisfies appropriate smoothness con- 
ditions, and that the statistics of the random delay 5 satisfy the relation 
that E[ | 8 — 8 \ k ] ^ h k X k ; where h is some constant, E denotes expecta- 
tion, 8 = E8 and X 2 is the variance of 8. The requirement on the random 
delay will always be satisfied when 8 is restricted to a bounded interval. 
Our model will also assume that intersymbol interference is negligible, 
or equivalent ly, that we are dealing with a single transmission. 

Under these assumptions an approximation is obtained for the test 
statistic which is optimum, in the Neyman-Pearson sense, for the case 
of gaussian receiver noise with exponential covariance. Generally the 
test statistic involves a nonlinear operation. However, for the case of 
phase reversal signals, only linear operations are required. 

The form that the test statistic takes for "white noise," which is con- 
sidered as a limit of the exponential covariance case, is obtained. It is 
shown that for phase reversal signaling the optimum receiver is a cross- 
correlator and that a portion of the correlation waveform is the expected 
value of the received signal itself. 

Fig. 1 depicts the communication system under consideration. A sig- 
nal, s (l) (t) (i = 1 or 2), which is non-zero only over an interval [0,7*1, 
is transmitted through a channel. The channel consists of a random delay 
8 and a linear time invariant filter whose output, x~ x) (t — 8), is disturbed 
by an additive noise source, n (t). It is assumed that the variance of the 
random delay, denoted by X 2 , is small. Furthermore, the noise, n(t), will 
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Fig. 1 — Model of Binary Communications System with a Small Random Delay. 
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be assumed to be a sample function of a stationary gaussian random proc- 
ess with mean zero and covariance R(t) = R (0)e~ p ' T '. For this system 
we seek the test statistic, based on an observation of the received signal 
over a fixed interval of length equal to the duration of the transmitted 
signal, which gives rise to the minimum error probability in the receiver 
decision process. 

II. DERIVATION OF THE TEST STATISTIC 

It is known that in testing between two simple hypotheses a Neyman- 
Pearson test will give rise to minimum error probability. Furthermore, 
Grenander 4 has shown that in the "regular case" the desired test 
statistic, which is a random variable called the likelihood function, I, 
can be obtained as the limit of an N dimensional likelihood ratio. In the 
subsequent development, in which it is assumed that we deal only with 
the regular case, the receiver test statistic is obtained as the limit of 
such an JV-dimensional likelihood ratio. 

The receiver input, y(t), is given by 

y(0-* w (l-«)+n(0, 1=1,2 (1) 

where x (i) (t — 5 ) is the portion of the input resulting from sending the 
signal s (i) (t) when the random delay is 5, and n(t) is the receiver noise. 
The noise is assumed to be gaussian with covariance R (t) = R (0)e~^ ' T '. 
Using a theorem due to Belayev 5 the noise sample functions can be 
shown to be almost surely continuous. Furthermore, almost all sample 
functions can be expanded almost surely in a pointwise convergent series 
in terms of the eigenf unctions of the noise covariance kernel. That is, 

n(t) = 2>*wb(0, (2) 

k 

with 

E(n/n k ) = <T?8 jk , 
where the <pk satisfy the integral equation 

<r*V*(0 = f ' du<p k (ti)R(t - u), to <t <h (3) 

J t 



and 



n k = I dt<p k (t)n(t). 
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Here t and h = t + T mark the beginning and end of the receiver proc- 
essing interval, 22 is used to denote J3 j the symbol E signifies mathe- 

k Jt=l 

matical expectation, and 

fl j = k 
8 jh = 

Since the {n*} are themselves gaussian and uncorrelated they are 
statistically independent. Assuming that the noise has zero mean, the 
joint density function of the first N coefficients, p N , can be written as 

P „(n il fc = l,..,i V )4n(i) ! iexp{-g ) (4) 

where the n* are ordered corresponding to the relationship that 

0"1 = 02 = ■ ■ ■ . 

Now consider a formal expansion of the receiver input in terms of the 
eigenfunctions of (3 ) . One can write 

y(t) ~ Z ym(t) = Z **%»*(<) + Z ****(*) 



k 



y k 4 ] dt<p k (t)y(t) = J dt<p k (t)[x s w (t) +n(t)] ( 5 ) 

** W) (S) A f dt<p k (t)x li) (t - 5), 

where by definition I = [t , ti]. Since the series Z nk<Pk(t) is almost surely 

pointwise convergent to n(t) we need only investigate the sense in which 
Z^* ( 'V*(0 converges to rc (,) (£ — 5). With this in mind we digress to 

k 

consider some of the properties associated with the eigenfunctions of the 
integral (3) with R(t — u) = R(Q)e~ p ' *~* '. It is easy to show that in 
this case the solutions of the integral equation are identical to those 
which satisfy the following differential equations and boundary condi- 
tions : 

d 2 ,, /3(/W - 2fl(0)) ,., _ 
ar erf 

(}<p k {U) = j t <p k (t)\ t ^ t0 (6) 

fiVkiti) = —j<p k (t)\i =tl . 
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The solutions of this system are proportional to 

cos y t [t - (*+*)] 

for k even, and 

sin y k [t - P±^)j 

for k odd. Here y k satisfies the relation (/3 2 + y k ) = 2/Stf (0 )/<r k . 

The differential equation (6) and the associated boundary conditions 
together form a Sturm-Liouville eigenvalue problem. The convergence 
properties of expansions in terms of the resulting eigenf unctions, the 
{<Pk (0 } , are stronger than those generally associated with expansions in 
terms of the eigenf unctions of the integral equation (3). An expansion 
of an integrable function f(t) on the interval (to , U) in terms of the 
eigenfunctions of a Sturm-Liouville system, possesses the following 
property: 6 

In every interval where f(t) is continuous and of bounded variation, 
the expansion converges uniformly and absolutely to f(t). If at the ends 
of the interval there are neighborhoods in which/ (t) is of bounded varia- 
tion then the series converges at these points to /(^o+) and f(h-). 

It will be assumed that the transmitted signal s (,) (t) and h(t) are such 
that rc (,) (t — 5) is continuous and of bounded variation. In fact to make 
the succeeding development valid we shall have to impose more stringent 
requirements on s (,) (t) and h(t). Under this assumption the expansion of 
x U) (t — 5) in terms of eigenfunctions of the integral equation (3) will 
converge uniformly and absolutely to .t (,) (t — 8). 

Returning to (5) we have established that ^2 ympkit) converges point- 
wise to y(t) for almost all sample functions. That is y(t) = 2 2/w*(0> 

k 

almost surely. 

Choosing the values of the {y k } set as the observable coordinates, the 
likelihood function, I, can now be determined as the limit of the ratio 
of two AT-dimensional density functions evaluated at the sample values, 
the \jj ,j = 1, • • • N. Here we have used the same symbol, y, , to repre- 
sent the sample and the random variable itself. 

Thus I can be written as 

l,h g^M , (7) 

»— v m (yi, •••,2/*) 

where p (,) (yi , • • • ,ys) is the joint probability density function of the 
first N members of the {y k } set. 

Noting that y k = ^a- (,) + n k and using the fact that the signal and 
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noise components are statistically independent the joint probability 
density of the first N of the y k can be written as a convolution. Thus 



P 



* m <Vi,~',Vw) = [ •• \ dz u ---,dz N S {i) (z u ---,z N ) 

J — CO •'— CO 



(8) 



■p N (yk - z k ,k = 1, • • -, N) 

pN U) (yi, ■•■,y N ) = E+(i)<pAy t - ** w \ k = 1,---,N) 

where 

J" (,) (21 , • • ■ , z N ) is the joint probability density function of the first 

N of the ^* (i) («), 
E+ (i) denotes an expectation with respect to the random vector, 

* W = (*i (i> («) ■■•*r W («)). 

Since the ^a- (,) (5) are all functions of the random variable 5, the aver- 
aging process can be performed with respect to 8 instead of the fa (8). 
Thus 

Pw W (jfi , • • • , Vw) - tf«Mlfc " fc W («), * - 1, • •' , N). (9) 
Using (4 ) , (7 ) and (9 ) , and cancelling common factors the likelihood 
function, I, is given by 



I = lim 

AT-»W 



E t exp I J2 


ZA- - 


2 




E t exp { 2 

U=i 


2/* - 


^■ (2, (5)" 
2 





(10) 



Let^^g^-^]^. 



Then (10) can be rewritten as 



I = lim 



. E s exp {Q N (1) [8]} 



(11) 



^coJ^exp {Q* (2 >[5]} ' 

At this point one may develop an integral form for Qh [8]. The in- 
tegral will suggest the form of lim Q N M [8] which will be required in the 

subsequent development. Define 



fc=l 0"* 

-V 



A-=l 

A' 

|fjr<0 = E^(0- 



(12) 
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Then using the orthonormality of the {*>*(0} one has the relation 

Q N (i) [8] = I dt [y N (t) - SM] q N (i) (f,S) . (13) 

It also follows from (3), (12), and the orthonormality relations, that 
q N (i) (t;8) satisfies the integral equation 

x „ (i) (t-8) = f du q N {i) (u;8)R(0) exp (-0\t- u |). (14) 

The solution of the integral equation is 

»"<*> -asAw^-S)* *- (15) 

Thus (13) can be rewritten as 

e.«w - / * [«.« - 2*M] ^(o) (" ! - Q *«*>. (i6) 

Let us return to the study of the likelihood function of (11 ). In general 
it seems that the expectations appearing in (11) cannot be evaluated. 
However it is possible to evaluate the required expectations when the 
variance of the random delay 5 is very small. With this in mind we drop 
the superscript and expand Q N (8) in a power series with remainder 
around 5 = 5. 

QM = <M*1 + (* - flOr'HI + (8 o. 5)2 Or' 1 * 1 

(17) 
+ (g 5) (?/"[5 + f(J - 5)], ^g = 1 

where we have used the notation 

Q*'® =^Q*(«)ki- 

It follows that exp {Q N [8]\ can be expressed as, 
exp {Q N [8]\ = exp {Q N [8]\ exp (5-5) a* 



( 3 - 5) 2 (5 - 5) 3 ,, 

+ 2 , &» + 3 , C "W 

where 

a* = Q/[5] 6* = Q/[8] c N (8) = Q„'"[5 + 0(5 - I)]. 
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On expanding the exponential in a powers series and averaging the 
uniformly convergent series term by term with respect to 5 one obtains 

E, exp { Q N [S] } - exp { Q N [I] } j"l + g *s] , (18) 

where 

d kN = E s [[(5 - 8)a N + ^A 2 b„ + ^-=J^ cAS) J 

The limit of the likelihood ratio of (10), Z, can be expressed as 

l = lim exp {Q„"> ® - Q„ m [a] | } '=H^Jr ■ (19) 

(»+s¥) 

Note that if for some number -D(JV) 

|4ur| < [#(#)]* for all fc, 

then the infinite sums each converge since the series is majorized by 
exp [D(N)\. It is observed that 

I 3i + 22 + 23 | fc ^ [3 max | Zf | ] fc 
which implies 

| ft + 22 + 23 |* ^ 3*( I 2l |* + I 2 2 |» + I 2 3 |*). 

Using the definition of d kN and the above inequality in conjunction with 
the Schwarz inequality yields 

I *» I £ 3* j| a N \*E | (5 - S) h | + |^|" # | (5 - 8) 2k | 

(20) 



+ u 



,c*(S) 



3! 



•^ U - 5 



i r) } • 



Now we restrict our investigation to the class of random delays whose 
probability distributions satisfy the relationship that for some number 
h and all k, 

A I («- 5)* | ^ fc*X*, (21) 

where X 2 is the variance of the distribution. This condition will be satis- 
fied by all probability distribution functions which take on the values 
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zero and one in a bounded region of the real line. That is when the values 
of 5 are essentially restricted to a hounded region. 
For distributions which satisfy (21 ) one finds 



|4ur| £8*{|*r|*A»X* + 



2! 



h X + 



j 3fc, 3* "j 

^- v i -an n*} ■ 



However if Q»" [8] is bounded, E \c N (8) \ 2k . ^ {c N ) 2 and 

I dun I ^ 3* { I Q„ | AX + | b N | fcV + (e*)W}*, I <fe* I ^ Am* 

where 



8* = sup Qn"[S], 
A H 4 3 { | a H | AX + | b N | fcV + (8*)W} . 

Now it can be shown that under certain restrictions on the channel 
and the class of transmitted signals the following limits exist almost 
surely : 

lim a N = a, 

N- » 



lim b N = b, 

N-oo 

lim c N (5) = c(8) 
limbic* (fi) T ^ c k 



(22) 



for all k. 

00 

The convergence of lim ^2 (d*jv/fc!) can now be demonstrated. In 

#-►00 fc=l 

this and the subsequent development we will consider convergence in 
the almost sure sense. Breaking the sum into a finite sum from k equal 1 
through m and a sum from m + 1 to °o and taking magnitudes gives 



oo , m 7 

Inn 2^ TT ~ 2-, n 



, . ^ dkN 

lim 2, -T7 
A '-oo fc-m+1 /C! 



where rft = lim dk# . Since | d k s | is less than or equal to some number 

A 1 -* 

An t then 



00 7 m I 

1,1,1 ^TT " ^- M 



AT-oo fc=,„+l A ! 



But lim A N = A exists, thus for all N > N 

AT-oo 
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\A - A N \ <t, | A H J < | A f + «, 



AT-.* fc-i /c! ^i fc! 






fc! 



By choosing m sufficiently large the right hand side of the inequality 
can be made less than any positive constant. Therefore 



For small X one makes the approximation 



lin. ("l + ± %fl 



1 + 4 + 



2!' 



(23) 



(24) 



where terms involving X to powers greater than X 2 have been neglected, 
and 

d k = lim dkN . 

Using the definition of dkN , which follows (18), we find 

di ~ (X 2 /2!)6, 

di ~ X"o . 

Restoring the superscript notation and using (22) and (24), (19) 
becomes 

Z~exp [limits] - Q N i2) [8]}\ 



•{1 + (X-"/2)[(a u, r - (ary + b w - b w ]\ (25) 

Equation (25) is the desired approximation for the likelihood func- 
tion for the case of small X. 



III. STRUCTURE OF THE OPTIMUM RECEIVER 

The structure of the approximation to the optimum receiver statistic, 
the likelihood function, can be determined from (25). The quantities 
appearing there can all be expressed in terms of the received signal, y(t) 
and the noise-free filter outputs x w (t) and x {2) (t), which are assumed to 
be known. In the Appendix it is shown that almost surely 
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■tfxit - 5) - x"(l-8)] dt 

- {/(In ~ 8) - 0x(/<, - 5)] L(/ ) - ^ 

+ [x(ti - 8) +(3x(tt - 8)] 



-5) 



(26) 



?y«.) - 



.t(/, - 6) 



■ [-ffx'it - 8) + .<•'"(/- a)] A 

- [.»•"(/, - 5) + fr'fa - «)] Mi) - ir(/, - *)] (27) 
+ [.r"(/ - I) - 13/ (t, - «)] |y(fc) - \X{U ~ I)] 
+ |[aj'(«, - 8) + 0x0, - 5)] [./(/, - 6)] 

- \[x{k - I) - Px(t - 8)] [z'0 o - 5)]|, 

■foV'0 - 8) - .r""(/ - 5)] rf* 

+ [0.r"(/„ - 8) - /"(/o - 8)] [//(/„) - .r(/o - 6)] 

+ \0.r"(ti - 8) + *'"(«, - 8)] [;/(/,) - xft - 5)1 (28) 

- fdt[(0x'(t - 8)) 2 - *'0 - ©*"'(« - 8)] 

+ fee"0o- 8) - /3.r'(/ - 5)].t'(/ - 8) 

- [.«•"(/, - 8) + 0z'0. - 8)]/0, - 5)1. 

It is also found that c(8) is a piece wise continuous function of 8. 

In obtaining these results it is assumed that the second derivative of 
x(t) is continuous and that the quantity [— 0V0 — 5) + x" (t — 8)] 
appearing in the integral in (27 ) is of bounded variation and continuous 
except at a finite number of points. Similar assumptions are made on 
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\f? x "(t - 8) - x"" (t - 5)] in the expression for b and a similar term in 
the expression for c(8). One requires that the fifth derivative of x(t) 
be of bounded variation and be continuous except at a finite number of 
points. These assumptions allow Sturm-Liouville expansions of the 
integrands involved in (27) and (28) which, using the orthonormality 
of the {<p k (t)\ , yield the proof of (27) and (28). Thus a sufficient 
condition for the validity of these results is that the second derivative of 
x(t) be continuous and that the fifth derivative of x(t) be of bounded 
variation and continuous except at a finite number of points. 

The desired test statistic, I, is obtained by substituting (26), (27) and 
(28) in (25) with x (i) (t), i = 1 or 2, replacing x(t) in the appropriate 
places. The resulting expression is nonlinear in y{t) and rather lengthy 
and will not be explicitly stated here. Observe however from (27) that 
if phase reversal signaling is used, that is x m (t) = —x w (t) , the quan- 
tity (a (1) ) 2 — (a (2) ) 2 will be linear in y(t). Furthermore under this condi- 
tion, to first order in X 2 , In I will be linear in y(t). Therefore for phase 
reversal signaling the receiver correlates y{t) with a signal related to 
v(t — 8) and its derivatives evaluated at 8 = 5. 

The "white noise case" will now be obtained as a limit of the exponen- 
tial covariance case. Setting R (0) = N f3/2 and letting -» w in the 
expressions for Q, a and b one obtains 



x(t - 8) 



+ ^r(fdt[y(t) - x(t - «)] x"{t - 8) 



(29) 



-f z dtW(t - 6)]*\\ 



:=*»> 



In the above expression the notation { j x=*(»> is used to indicate 
that the expression in braces is evaluated with x(t) = x {2 (t) and the 
result is subtracted from the result which obtains when x(t) = x ll) (t). 
For phase reversal signalling, x w (t) - -x (i) (t) = x(t), (29) becomes 



- x\U - «)) x'(l -8)- x"( - 5)11 . 



(30) 
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Let us examine the correlation waveform appearing inside the braces 
in (30). The first term represents the receiver input when the random 
delay is equal to its mean value, 8 = 5. If X is set equal to zero the re- 
maining terms vanish and one obtains the standard result which is that 
the receiver input should be correlated with the signal portion of the 
input, x(t — 8). The remaining terms inside the braces in (30) are the 
perturbations introduced by the random delay. 

Let us expand x (t — 8) in the Taylor series with remainder 

x{i - 8) = x(t - 8) + (8 - 8) [-x(t - 8)] 

+ [(5 - 5) 2 /2!] [x"(t - 8)] + (Jt(S - 8,0, 

If x'"(t — 5) is continuous and the moments of 8 satisfy 

E s | (5 - 8) k | ^ h k \ k , 

then for small X, x(t — 8) + (K 2 /2)x"(t - 5) is the principal part of 
Esx(t — 8). Thus, part of the correlation waveform is essentially the ex- 
pected value of the received signal. The other term in the correlation 
waveform involves x (t — 8), which is the slope of the received signal 
when the delay is 5. The weight attached to it is proportional to the 
difference in the squared values of the received signal at time t and at 
time h when 8 = 8. As yet no physical significance has been found for 
this term. 

IV. CONCLUSIONS 

An approximation has been obtained for a test statistic that minimizes 
the error probability of a binary communication system which operates 
over a linear channel, with a small random delay, in the presence of gauss- 
ian noise of covariance R(0)e~ p ' T '. For the case of phase reversal sig- 
naling, the statistic, which in general is nonlinear, is a linear functional 
of the receiver input. Treating the "white noise" case as a limit of the 
exponential covariance case, the test statistic is expressed as a cross-corre- 
lation operation. The waveform with which the input is correlated is 
related to the expected value of the received signal plus a term propor- 
tional to the slope of the received signal when the delay is equal to its 
mean value. 

V. ACKNOWLEDGMENT 

The author wishes to express appreciation to Dr. A. R. Cohen of 
New York University, Dr. I. Jacobs of Bell Telephone Laboratories and 
the referees for their valuable suggestions. 



1136 THE BELL SYSTEM TECHNICAL JOURNAL, JULY-AUG. 1965 
APPENDIX 

The Evaluation of Terms Arising in the Approximation to the Likelihood 
Function 

The convergence asserted in (22) will now be established and the 
structure of the terms appearing in the likelihood function, (25), will be 
exhibited. 

First let us note the following property associated with an expansion 
in terms of the eigenfunctions of a Sturm-Liouville system. 

Let g(t) be a piecewise continuous function which is of bounded varia- 
tion and let z(t) be a piecewise continuous function. If the sets {gh\ and 
{zk\ are the expansion coefficients of g(t) and z(t) in terms of the above 
eigenfunctions, then 

[ dtg(t) z(t) = Eflrffc. (31) 

This expression is obtained by noting that the series expansion of g(t) 
converges uniformly except at a finite number of points. 

One considers the following integral which is suggested by (16), 

*«> - mm l dt [*<* ~ t] if - S) •<* - o. ™ 

where one requires that [/3 2 — (3 /dt )]x(t — 5) be of bounded variation 
and piecewise continuous. Let 

g{t - 8) 4 (? - Q xit - S) = fait - 8) - x"it - 8}. (33) 

git) can be expanded in a series in terms of the eigenfunctions of (6). 
This series converges uniformly to git) where git) is continuous and to 
\[giu + ) + fif(w-)] at points where git) is not continuous. Thus 

00 

git -8) = YsQWkU), (34) 

4=1 

where 

g k = J du ip k iu) git — 8). 

By integration by parts and utilizing (6) one can show that 

9k = iP 2 + 7k )h + <Pk(U)[x (U - 8) - fix'iU - 5)] 

(35) 
-<Pkiti)[x ih - 8) + &x{h - 8)), 
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where 

\pk is the fcth expansion coefficient of x(t — 8) in terms of the \<p k (t)} 
set, 

tf + 7*') = pPB(0)/n*\, 

<Tk is the fcth eigenf unction of (6). 

Substituting the expansion for g(t) in (32) and using (35) one finds 

-I- [x'(t - 8) - Px(U ~ 8)] g (y* - ^J *>*(*>) (36) 

- [x'(h - 8) + px(k - «)] g (y* " ^)f*(ft)}. 

Noting that X ( 2A — 77 ) — , equals lim Qa-[5] one finds 

r n m jr/sN W (to - 8) - 0x (to - 8)] ^ / ^A , . , . 

hm Q N [8] = K (8) 2BR(0) ^ V* ~~2/ <Plc ™' ^ ' 

, [x ( h- 8) +0x(h- 8)} A / *A . 

+ mm S v" " 2) vM 

One now can proceed to evaluate a — lim qn . 

AT=0O 

From the definitions following (17) and (10) one has 

a N = QA8] = t (y »-?* (g)) x k ' (8) . (38) 

i = l Ok 

Consider 
K:(5)A J[ <fc [y(0 - X(t ~ 8) ~\ l-P 2 *(t ~ V + *"'« " «)]• (39) 

Let us assume that [—fx(t — 8) + a/"(< — 5)] is piece wise continu- 
ous, [fix (I — 8) — x" (t — 8)\ is continuous and that both are of bounded 
variation. Define 

a* = J dt <p k (t) [-{?x(t - 8) + x'"(t - «)]. (40) 

Then under these assumptions and using (33) and (34) one finds 

a k = g k (8). (41) 
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Applying (31) to Ki(8) gives 



*i(«) -g(»-^)*- (42) 



Using (41), (42), (35) and (31) yields 

+ m ^ - S) - ^* - « ] [» W " ^] (43) 
- \ [x(to - 8) - fix (U, - S)]x(to - 8) 

+ § fr'tti - 5) + /3a: (h - 8)]x(h - 8). 
The left hand side of (43 ) evaluated at 5 = 5 is lim a N . 

JV-» oo 

Proceeding in a similar manner the lim b N and lim c N (8 ) are ob- 

AT-»oo JV-»oo 

tained. 
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